First-principles analysis of spin-disorder resistivity of Fe and Ni 
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Spin-disorder resistivity of Fe and Ni and its temperature dependence are analyzed using non- 
coUinear density functional calculations within the supercell method. Different models of ther- 
mal spin disorder are considered, including the mean-field approximation and the nearest-neighbor 
Heisenberg model. Spin-disorder resistivity is found to depend weakly on magnetic short-range or- 
der. If the local moments are kept frozen at their zero-temperature values, very good agreement 
with experiment is obtained for Fe, but for Ni the resistivity at elevated temperatures is signif- 
icantly overestimated. Agreement with experiment for Fe is improved if the local moments are 
iterated to self-consistency. The overestimation of the resistivity for paramagnetic Ni is attributed 
to the reduction of the local moments down to 0.35^_b. Overall, the results suggest that low-energy 
spin fluctuations in Fe and Ni are better viewed as classical rotations of local moments rather than 
quantized spin fluctuations that would require an [S S correction. 



I. INTRODUCTION 

Electron scattering off of spin fluctuations in magnetic 
metals results in an "anomalous" contribution to electric 
resistivity. ^'^''^ The analysis of this spin-disorder resistiv- 
ity (SDR) is of interest because it can provide material- 
specific information on the character of spin fluctuations 
which is not easily accessible by other means. Scattering 
on spin disorder is also an important factor degrading the 
performance of magnetoresistive nanostructurcs in spin- 
tronic devices. 

The magnitude of the spin-disorder contribution to re- 
sistivity is comparable to the phonon contribution near 
and above the Curie temperature Tc-^ (Magnetic scatter- 
ing ampHtudes have no small parameter unless the local 
moments are small.) It is usually assumed that SDR 
is constant well above T^- In this region Matthiessen's 
rule is valid, and the phonon contribution can be fitted 
to the Bloch-Griineisen formula. The excess resistivity 
in the whole temperature range may be attributed to 
spin disorder)^ although one may expect deviations from 
Matthiessen's rule at low temperatures where transport 
is carried by weakly interacting spin channelsi^ In addi- 
tion, it was argued that in some cases (such as Ni) spin 
disorder may change the character of states on the Fermi 
level and thereby appreciably change the phonon contri- 
bution itself.^'- 

Many authors have studied SDR theoretically using 
the s-d model Hamiltoniani^iii^i^ In this approach the id 
shells in transition metals (or / shells in rare earth metal- 
lic magnets) are assumed to be localized at atomic sites 
and partially filled, forming magnetic moments that 
are coupled to the current-carrying conduction electrons 
by exchange interaction Hsd = ~JsdYlii SiS^, where Jsd 
is the local s-d exchange coupling constant and the 
spin-density operator of the conduction electrons at site 



i. Thermal fluctuations of the d-electron spins generate 
an inhomogeneous exchange potential; in the Born ap- 
proximation the SDR is then determined by the conduc- 
tion electron band structure and the spin-spin correlation 
functions of d-electron spins If the scattering is approx- 
imated as being elastic, only equal-time spin correlators 
have to be considered. Further, if the mean-field approx- 
imation (MFA) is used for 2>d spin statistics, the SDR 
behaves as Pmag(r) = po[l - NP{T)/S{S + 1)], where 
M{T) = (S(T)) is the magnetization at temperature T 
and po (X J'^j,S{S + 1).^ Note that above SDR is con- 
stant and equal to po- The shape of the Fermi surface of 
conduction electrons is immaterial to this prediction as 
long as the scattering is elastic 

The effects of magnetic short-range order (MSRO) 
on SDR have also been investigated within the s-d 
model-i^iiSiiiii^ii^ii^ii^ Thls problcm has attracted con- 
siderable attention in connection with a "bump" in the 
resistivity that is observed near in some magnetic 
metals (although it is usually quite small). ^ The anal- 
ysis of critical MSRO effects showed that a cusp may 
appear near Tc due to long-wave critical fluctuationSfi 
although it should usually be strongly suppressed by fi- 
nite mean-free path and cancelations due to Fermi surface 
integration, y. It was also found that the effect of MSRO 
and even its sign are sensitive to such details of the model 
as the conduction band occupation and the form of the 
scattering (pseudo)potential .^^'^'^'^'^ 

Although the s-d model provided physical insight into 
the mechanism of SDR, it suffers from serious limitations. 
First, the distinction between localized and conduction 
electrons is not justified in transition metals where Zd 
electrons are itinerant and form the Fermi surface. Even 
if the current is dominated by light s-like bands that 
can be distinguished from heavy c?-like bands, the relax- 
ation rate is dominated by interband {s-d) scatteringii^ 
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Second, at elevated temperatures the scattering potential 
generated by spin disorder is of the order of the exchange 
splitting, which is not small compared to the bandwidth. 
This invalidates the Born approximation which is usually 
made in model calculations. Third, the s-d model does 
not properly take into account the change of electronic 
structure due to disorder. 

The first-principles approach to SDR is free from all 
these limitations and can be used for quantitative calcu- 
lations of SDR. This is of particular interest for the the- 
ory of itinerant magnets, because, as mentioned above, 
SDR depends on spin-spin correlation functions. Differ- 
ent theories of itinerant magnetism make conflicting pre- 
dictions for such properties as the degree of MSRO, the 
mean-squared magnetic moment, and their temperature 
dependenceii^ii^i^i^Si^ these quantities are quite hard 
to measure directly. By calculating SDR for a particu- 
lar model of spin fluctuations and comparing the results 
with experiment, one can attempt to validate or rule out 
different spin- fluctuation models. 

Earlier we have calculated the temperature dependence 
of SDR in Fe and Ni using supercell calculations within 
the tight-binding linear muffin-tin orbital (TB-LMTO) 
method using the mean-field distribution for spin orien- 
tation statistics. Good agreement with experiment was 
obtained for Fe, but for paramagnetic Ni the SDR was 
found to be significantly overestimated. In this paper 
we analyze the temperature dependence of SDR for Fe 
and Ni in greater detail and consider the effects of mag- 
netic ordering, MSRO, and local moment reduction. We 
also study the importance of the basis set size and self- 
consistency of the atomic potentials. 



II. GENERAL APPROACH AND METHODS 

Our approach is based on noncoUinear density func- 
tional theory (DFT). All the valence electrons are treated 
on the same footing, while the scattering potentials are 
determined by the self-consistent electron charge and 
spin densities. We use the TB-LMTO method^^ which 
represents the electronic density of the crystal as a su- 
perposition of overlapping atomic spheres; the electronic 
density inside each sphere is spherically symmetric. This 
method is known to work very well in close-packed mate- 
rials, and it allows us to introduce spin disorder in various 
ways. In this work we used the rigid spin approxima- 
tion which assumes that the spin density in each atomic 
sphere remains collinear, while the spin densities of dif- 
ferent atomic spheres become noncollinear at finite tem- 
peratures. In the simplest model the electron charge and 
spin densities in all atomic spheres are taken from the 
ground state and frozen, while the directions of the spin 
moments in different spheres are randomized with the an- 
gular distribution function taken from MFA at the given 
temperature. This model is expected to work reasonably 
well for Fe which has a fairly stable local moment In 
Section mil we show that this is indeed the case; however, 



for Ni the paramagnetic SDR calculated in this way is 
about twice too large. In order to explain this discrep- 
ancy, the dependence of SDR on the degree of MSRO 
and on the magnitude of the local moment is studied in 
Sections HV] and El 

We use the supercell approach and calculate the areal 
conductance of a layer of spin-disordered metal FM(D) 
sandwiched between fully ordered semi-infinite leads 
FM(0) made of the same metal (see Fig.[T]). The resistiv- 
ity is then proportional to the slope of the inverse conduc- 
tance as a function of the disordered layer thickness, once 
the Ohmic limit is reached. For the given thickness of 
the FM(D) layer, the conductance of the system was av- 
eraged over several disorder configurations (typically 15). 
The planar system is represented by a laterally periodic 
prism with an axis along the [001] crystallographic direc- 
tion, and care is taken to make sure that the conductance 
scales as the cross-section of the prism. To calculate the 
conductance we use the principal-layer Green's function 
techniqu o^^'^^ and the Landauer-Biittiker formalismSl in 
the implementation allowing for noncoUinearity in the 
active region;^ This technique was employed before to 
study the effects of substitutional disorder on transport 
in magnetic multilayers;^ it is similar to the supercell 
Kubo-Greenwood method used to calculate the residual 
resistivity of binary alloysi^ 



L lead 


active region 


mmA 


T T 
T t 
T 1 


T t T 

T r t 

T T T 


J \ ^^-*\^ 


t T T 
T T T 
T T T 


T r 
t r 









FIG. 1: (Color online) The schematic picture of the system 
used in the calculations. Vertical lines indicate the embedding 
planes. 



If the atomic potentials in the supercell are not con- 
verged to self-consistency with the given spin disorder 
configuration, care needs to be taken to ensure local 
charge neutrality. Indeed, FM(D) and FM(0) materi- 
als have different Fermi levels that must normally be 
matched by the contact voltage. In order to enforce 
charge neutrality in the FM(D) region, a constant po- 
tential shift was introduced in this region so that the 
charge in the central part of FM(D) averaged over disor- 
der realizations was zero. This potential shift plays the 
role of the contact voltage. Note that no matter how the 
FM(0)/FM(D) interfaces are treated (self-consistently or 
not), they add contact resistances to the circuit. How- 
ever, since the resistivity of the FM(D) material is ex- 
tracted from the thickness dependence of the resistance 
in the Ohmic limit, the simplified treatment of interfaces 
has no effect on the results. 
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III. SPIN-DISORDER RESISTIVITY IN THE 
MEAN-FIELD APPROXIMATION 

A. Paramagnetic state 

In this section we analyze the temperature dependence 
of SDR for iron and nickel using MFA for thermal spin 
disorder; the spin-spin correlator is purely local in this 
approximation. First we consider the paramagnetic state 
where the angular distribution function is isotropic, and 
the resulting SDR is temperature- independent. 

We need to make a physically reasonable choice of 
atomic potentials for the conductance calculations. It is 
known that the local moments in Fe are quite stable;^^ in 
particular, the DLM method, which employs the coher- 
ent potential approximation for spin-disordered states, 
shows only a small reduction of the local moment in para- 
magnetic Fe compared to its ground-state valuci^ As 
seen below, direct averaging of self-consistent local mo- 
ments in the paramagnetic states gives a similar result. 
Therefore, for Fe it is reasonable to use frozen atomic 
potentials taken from the zero-temperature ground state 
in all calculations. We have also checked the effect of 
self-consistency on SDR in Fe and found it to be small 
(see below). The situation is entirely different for Ni, 
where the local moment depends on the magnetic state; 
in particular, it vanishes altogether in the paramagnetic 
DLM approximation.'*'* Since longitudinal spin fluctua- 
tions (that are absent in our approach) can at least par- 
tially restore the local momentsj^ it is not a priori clear 
how the atomic potentials should be modified for Ni. In 
this section we use frozen atomic potentials; the neces- 
sary corrections are discussed later. 

Fig. [2] shows the inverse areal conductance for param- 
agnetic Fe and Ni as a function of the disordered FM(D) 
region thickness. Here we used the frozen ground-state 
atomic potentials and the LMTO basis including s, p, 
and d orbitals [Imax — 2). The supercell cross-sections 
contained 4x4 (for Fe) and 3x3 (for Ni) cubic unit cells 
with edges oriented along the [100] directions. Almost 
perfect Ohmic behavior is apparent for both Fe and Ni, 
which establishes the validity of our approach. 

Table [J lists the values of SDR found for paramagnetic 
Fe and Ni using different supercell cross-sections, LMTO 
bases truncated at Imax = 2 and Imax — 3 (the latter 
includes / orbitals), as well as the value found using self- 
consistent (rather than frozen) atomic potentials for Fe. 
It is seen that the results are well converged with respect 
to the supercell cross-section, and even 2x2 supercells 
provide sufficient accuracy. This is reasonable because 
the mean-free path in the paramagnetic state is small. 

The calculations with self-consistent atomic potentials 
were performed as follows. In order to reduce the statisti- 
cal error, the averaging of the conductance was performed 
using the same sets of random spin disorder configura- 
tions as in the calculation with frozen potentials. For 
each individual spin configuration the atomic potentials 
were iterated to self-consistency using the Fermi distri- 
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FIG. 2: (Color online) The area-resistance product AR of the 
FM(0)/FM(D)/FM(0) systems as a function of the FM(D) 
layer thickness for bcc Fe (black filled symbols) and fee Ni (red 
or gray empty symbols) obtained with Imax = 2. Circles and 
squares correspond, respectively, to the paramagnetic state 
and to the lowest temperature for which the calculations were 
made. 4x4 and 3x3 supercells were used for Fe and Ni, 
respectively, with edges along the [100] directions. Straight 
lines show the linear fitting; error bars are smaller than the 
size of the symbols. 



TABLE I: Spin-disorder resistivity in /ifi-cm for paramagnetic 
bcc Fe and fee Ni. The calculated values are given for basis 
sets with Imax = 2 and 3, as well as for different lateral cell 
sizes with edges along the [100] directions. SC denotes calcu- 
lations with self-consistent potentials. Standard deviations of 
SDR due to limited disorder sampling are included. 
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bution function corresponding to the experimental Tc of 
Fe. The resulting distribution of the sites over the mag- 
nitude of the local magnetic moment is shown in Fig. [31 
this distribution is Gaussian with a rather small width. 
The average local moment is only reduced by 3-4% from 
its ground state value. This small reduction appears 
to be similar to the DLM calculations of Ref. ^201, while 
Ref. !24 obtained a somewhat larger reduction. The self- 
consistent density of states (not shown) is very similar 
to the one generated by the frozen ground-state atomic 
potentials (see Fig.^J; below). 

The addition of / orbitals to the LMTO basis reduces 
the calculated SDR by approximately 15% for both Fe 
and Ni. Self-consistency in the paramagnetic state of Fe 
results in a similar reduction. This similarity suggests 
that the main reason for this SDR decrease is the reduc- 
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FIG. 3: (Color online) Distribution of the local magnetic 
moment in self-consistent fully spin-disordered bcc Fe. The 
Fermi temperature is equal to the experimental Tc. The ver- 
tical line shows the local moment at T = 0. The red (solid) 
curve shows the Gaussian fit to the data. 

tion of the local moment, which is, incidentally, very sim- 
ilar in both cases. In order to check this, we performed 
additional calculations for Fe in which the / channel was 
added to the basis while the charge density was kept un- 
changed from the self-consistent one with Imax = 2. For 
the frozen potential case, SDR reduced slightly from 106 
to 100 /ifi-cm; for the self-consistent paramagnetic case, 
it only reduced from 88 to 86 /xfi-cm, which is within the 
error bar. Thus, the effect of Imax pe-r se on SDR is very 
small for Fe. This is somewhat different from the binary 
alloy systems considered by other authors using both TB- 
LMTO and KKR (Korringa-Kohn-Rostocker) methods, 
where a larger effect of adding / states was foundi^i^^ 
In view of the weak dependence of SDR on Imax, below 
we use Imax = 2 in all calculations for T < Tc- 

The experimental estimates of SDR in the paramag- 
netic state** are listed in the last column of Table HI The 
agreement with experiment for Fe is quite satisfactory, 
and it is in fact improved if the reduction of the local mo- 
ment is included. In Ni the SDR calculated with frozen 
atomic potentials is overestimated by a factor of 2. This 
is not surprising, because, as mentioned above, the use 
of frozen atomic potentials is not justified for Ni. In or- 
der to understand the origins of the disagreement with 
experiment for Ni, possible modifications of the statisti- 
cal model for the paramagnetic state must be considered; 
this is done below in Sections IV and V. 

Recently, Buruzs et al^ calculated the SDR for Fe and 
Co using the disordered local moment (DLM) approach 
within the Korringa-Kohn-Rostocker method and found 
that their method significantly overestimates the param- 
agnetic SDR in these metals. The source of disagreement 
with our supercell method for Fe is unknown to us. 



B. Temperature dependence in the ferromagnetic 
state 

In this section we consider ferromagnetic state of Fe 
and Ni. We use frozen ground-state potentials and the 
basis with Imax = 2. As mentioned above, this approxi- 
mation is reasonable for Fe, while for Ni it is not appli- 
cable at high temperatures; nevertheless, comparison of 
these two systems will allow us to draw important conclu- 
sions. For the ferromagnetic state the spin configurations 
were generated using the mean-field distribution function 

where 9 is the angle between the local moment /x and 
the magnetization axis, M[T) is the magnetization at 
temperature T in MFA, and Hoff is the Weiss field. This 
distribution function depends only on T /Tc- 

Before we turn to the temperature dependence of SDR, 
let us look at the electronic structure of Fe and Ni with 
spin disorder. The spin-resolved DOS of Fe and Ni is 
shown in Figs. S] and [5] for several temperatures. These 
data were obtained by projecting the site-resolved DOS 
onto local spin-up and spin-down states (in the local ref- 
erence frame where the z axis is parallel to the local 
moment) and subsequent averaging over bulk-like sites 
and spin disorder configurations generated according to 
Eq. ([1]). The paramagnetic DOS of Fe is very similar to 
the KKR-DLM results.^'* As the temperature is increased 
from to Tc, the spin-up and spin-down states ran- 
domly hybridize with each other, the peaks broaden, and 
the van Hove singularities are washed out. The mean- 
squared deviation of the DOS from its average (shown by 
dashed lines) is quite small, which is a direct consequence 
of the large coordination number. In Fe the spin splitting 
is almost independent on temperature, while in Ni it is 
much reduced as T gets close to Tc- Note that the frozen 
atomic potentials in Ni are very far from self-consistency 
at elevated temperatures, but a self-consistent treatment 
neglecting longitudinal spin fiuctuations would be mean- 
ingless. We will return to this issue in Section V. 

Let us now discuss the temperature dependence of 
SDR. While we found above that 2x2 supercells were suf- 
ficiently large for the paramagnetic state, additional care 
needs to be taken at lower temperatures where the mean- 
free path becomes longer. We found that 4x4 supercells 
for Fe and 3x3 for Ni demonstrate linear dependence 
of the conductance on the length of the supercell for all 
temperatures down to about Tc/i (see Fig. [5]). This be- 
havior agrees with a simple mean-free path estimate us- 
ing the free-electron formula I = |^i?bai/p, where Ai?bai 
is the ballistic area-resistance product; / does not exceed 
the lateral cell size in this temperature range. Another 
indication of the Ohmic behavior comes from the distri- 
bution of the current over the spin channels. The conduc- 
tance of the FM(0)/FM(D)/FM(0) system is a sum of 
four partial conductances, G||, Gxi? G^j, G^f (the latter 
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FIG. 4: Spin-resolved density of states (solid lines) for bcc 
Fe averaged over random spin configurations with the mean- 
field distribution function HJ; (a) T = 0, (b) T = 0.25Tc, (c) 
T = O.STc, (d) T = 0.75Tc, and (e) T = T^. Dashed lines 
show the mean-square deviation of the DOS on a given site 
from its ensemble average. 



two are equal). Spin-conserving and spin-flip scattering 
have similar rates in our spin-disorder problem (as long 
as the temperature is not too low), and therefore the 
electrons "forget" their spin over their mean-free path. 
Therefore, in the Ohmic limit the partial conductances 
must be proportional to the number of channels in the 
left and right leads for the corresponding spin channels: 
Gctct' oc M^M^. This implies that in this regime we 
should have G-^-^Gn = GuGip This relation does in- 
deed hold down to T ~ Tc/3 unless the thicknesses of 
the FM(D) region is very small. 

The dependence of the calculated SDR for Fe and Ni 
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FIG. 5: Same as in Fig. |4]but for fee Ni. 



on the magnetization is plotted in Fig. [6] along with the 
experimental data^- (those for M{T) were taken from Ref. 
1371 ). The results for Fe agree rather well with experiment, 
especially at lower temperatures where the magnetic ex- 
citations are dominated by spin waves and our classical 
approach is, strictly speaking, invalid. This surprising 
finding is due to the fact that SDR in Fig. [5] is plotted as 
a function of the long-range order parameter and that, 
as we show below in Section llVi the SDR in Fe is in- 
sensitive to MSRO. The calculated SDR exhibits linear 
dependence on Af^(T) up to Tc, while the experimen- 
tal data deviate downward from the straight line. This 
deviation may be attributed to a small reduction of the 
local moment at elevated temperatures, as discussed in 
the previous section. 

For Ni the deep low-temperature region could not be 
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FIG. 6: (Color online) Dependence of spin-disorder resistivity 
on the magnetization for (a) Fe, and (bj^ Ni. Black circles de- 
note experimental data combining Ref. 4 for pmag(r) and Ref. 
[37I for M{T). Blue (black) squares show mean- field calcula- 
tions, filled red (gray) triangles denote Monte Carlo results, 
and green (gray) diamonds show reverse Monte Carlo calcu- 
lations. The empty red (gray) triangles show Monte Carlo 
results with larger cells: 6x6 for Fe and 4x4 for Ni. The 
upper axis shows temperatures corresponding to the given 
magnetization in MFA. The error bars along the x and y axes 
show statistical uncertainties where they exceed the size of 
the symbols. All results are for Imax ~ 2. 



accessed due to the increased mean-free path. Still, 
the agreement with experiment at lower temperatures is 
good, while at higher temperatures the calculated SDR 
strongly deviates upwards from experimental data. This 
deviation indicates the inadequacy of our spin fluctua- 
tion model; its possible modifications are studied in the 
following sections. 

The qualitative features of the calculated temperature 
dependence of SDR (with frozen atomic potentials) are 
different for Fe and Ni. It is seen in Fig.[n]that for Fe the 
SDR is proportional to 1 - AP{T)/AP{0) in agreement 
with the predictions of the s-d model if spin fluctuations 
are treated classically. On the other hand, for Ni this 
relation does not hold. As mentioned in the Introduction, 
the change of electronic structure due to spin disorder 
may lead to deviations from s-d model predictions. 

As seen in Figs. [Hand O the densities of states change 
quite appreciably with temperature for both Fe and Ni. 
Therefore, it may seem surprising that for Fe the tem- 
perature dependence of SDR agrees with the s-d model. 
Still, one can understand the difference between Fe and 
Ni using the following considerations. First, as seen in 
Fig. [5l the exchange splitting in Ni is strongly reduced 
at elevated temperatures, which results in the lifting of 
the heavy majority-spin 3d bands up to the Fermi level. 
Scattering into these flnal states from the light bands 
becomes possible, which decreases the lifetime of the lat- 



ter. This mechanism was invoked by Mott^ to argue that 
the reduction of the spin splitting in Ni can result in an 
anomalous temperature dependence of the phonon resis- 
tivity. The same argument applies to SDR considered 
here. According to Fig. [S] this happens approximately 
at T = 0.75Tc, which roughly corresponds to the upturn 
of SDR seen in Fig.lHb. On the other hand, for Fe, as seen 
in Fig. m the exchange splitting is constant and no new 
bands appear at the Fermi level. Consequently, no ad- 
ditional temperature dependence is introduced and SDR 
scales as 1 - M^{T)/M^{0). 

While plausible, the above arguments are not conclu- 
sive, because they assume without proof that the scatter- 
ing matrix elements between the light and heavy bands 
are large. On a more subtle level, one may speculate 
that the difference between Fe and Ni can be under- 
stood based on the relation between disorder broadening 
and spin splitting. At the given wavevector, the spectral 
function consists of delta-function peaks corresponding 
to majority and minority-spin states. In the presence 
of spin disorder, the spin states on neighboring sites are 
allowed to hybridize with random matrix elements, and 
the delta-function peaks broaden. At low temperature 
the broadening is small, and the peaks corresponding to 
different spins are well separated in energy from each 
other. However, at higher temperatures some of these 
peaks can merge and form common, "virtual-crystal-like" 
bands. Calculations of the paramagnetic spectral func- 
tions using the DLM method indicate that in Fe the ma- 
jority and minority-spin states remain separated through 
large portions of the Fermi surface even above Tc,"^^ On 
the other hand, in Ni the majority and minority-spin 
states are mixed in the paramagnetic state. Therefore, 
at certain temperature below Tc there is a crossover from 
separated to mixed-spin bands. The lifetime is expected 
to decrease as the bands merge, which again explains the 
upturn of SDR from the straight line in Fig. 



IV. EFFECT OF MAGNETIC SHORT-RANGE 
ORDER 

As mentioned above, short-range order can sometimes 
have a significant effect on resistivity. In this section 
we analyze the effect of MSRO on SDR in Fe and Ni. 
In particular, it is important to check whether the large 
disagreement with experiment for Ni found in Section Hill 
can be due to the use of MFA which neglects MSRO. 
This is especially interesting because strong MSRO in 
Ni has been suggested by some experimentsii^ and 
theories ji^ii^ 

Spin disorder configurations with MSRO were gener- 
ated using the Monte Carlo (MC) method for the classical 
Heisenberg model with nearest-neighbor (NN) exchange 
interaction on bcc and fee lattices (for Fe and Ni, re- 
spectively). These configurations were used to calculate 
SDR as described above, which can then be compared 
with MFA results. As before, we usually used 4x4 and 
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TABLE II: Spin-spin correlators Coi for the first three shells of nearest neighbors (i = 1,2,3) and the local correlator (i = 0) 
are shown for Fe and Ni for each considered temperature in Monte Carlo simulations and for reverse Monte Carlo method. The 
values of the correlators were found using 4x4 cells and averaging over the lengths that were used in transport calculations. The 
Curie temperature for bcc and fee nearest-neighbor Heisenberg model was obtained using the fourth order cumulant method. 
The values of SDR are compared with MFA results corresponding to the same . The listed uncertainties are due to the 
hmited disorder sampling. 
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3x3 supercells for Fe and Ni, respectively, which was 
sufficient to achieve Ohmic scaling of the conductance. 

The results of SDR calculations for this model are 
shown in Fig. [HI The MC results for SDR are very close to 
MFA results corresponding to the same magnetization, in 
spite of the presence of MSRO in the MC model. We also 
performed additional calculations with larger cell cross- 
sections (6x6 for Fe and 4x4 for Ni) in order to check 
whether the finite-size effects are important close to Tc- 
The results are shown by empty and full inverted trian- 
gles in Fig. [HI the temperature was taken to be the same 
as for the neighboring point for a smaller cross-section. It 
is clear that finite-size effects have no appreciable effect 
on SDR. Thus, is appears that the spin-spin correlations 
on the length scales comparable to our supercell lateral 
dimension have negligible effect on SDR, in agreement 
with the conclusions of Ref. 10. 

These results clearly show that MSRO characteristic 
for the NN Heisenberg model has virtually no effect on 
SDR in Fe and Ni. The magnitude of MSRO is illus- 
trated in Table [IT] where the spin-spin correlators for the 
first three shells of nearest neighbors are shown together 
with the corresponding resistivities. While MSRO in the 
NN Heisenberg model for the close-packed lattices con- 
sidered here is not strong, it is seen that its effect on SDR 
is much smaller even compared with the values of the NN 
spin-spin correlators. This insensitivity is likely due to 
the averaging over all the electronic states on the Fermi 
surfacCfi^, which should be very effective in destroying 
the interference from scattering at different sites in tran- 
sition metals with complicated Fermi surfaces. In fact. 



this averaging is also responsible for the small standard 
deviation of the local DOS from its mean (Fig. [4]) and 
justifies the DLM approach for transition metals. 

The spin-spin correlation function in real materials 
may be more complicated than in the NN Heisenberg 
model. However, if the interaction has a longer range 
while remaining mainly ferromagnetic, the MSRO must 
be weaker compared to the NN modeli^ First-principles 
calculations for both ferromagnetic and paramagnetic 
nickel show that the exchange parameters beyond nearest 
neighbors, while being much smaller than the dominant 
NN exchange, stay mainly ferromagnetiC] ^°'^^ Interaction 
of this kind can not support stronger MSRO compared 
to the NN Heisenberg model. 

Since the complete insensitivity of the resistivity to 
MSRO within the NN Heisenberg model is somewhat 
surprising, we have checked whether it is possible to ob- 
serve some change of SDR using spin configurations with 
an artificially introduced stronger MSRO. For this pur- 
pose we used the reverse Monte Carlo (RMC) method^ 
to produce a set of spin configurations with almost zero 
magnetization and deliberately targeting strong MSRO 
in the NN shell. Due to geometrical constraints, the spin- 
spin correlators in different neighbor shells are not inde- 
pendent. We found it quite difficult to produce strongly 
correlated nearest neighbors and at the same time avoid 
unphysical artefacts in the long-range behavior of the cor- 
relation function. The spin-spin correlators for the first 
three shells of neighbors in our RMC model are listed in 
Table [ITl The corresponding values of SDR calculated 
for Fe and Ni with this set of spin configurations are also 
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listed in Table HIl and shown by full and empty diamonds 
in Fig. [S] Here we used 4x4 supercells for both Fe and 
Ni and checked for finite-size effects using 6x6 supercells 
for Fe (essentially no difference was observed compared to 
4x4 cells). As seen in Table III the MSRO in this model 
is significantly stronger compared to the NN Heisenberg 
model. The effect of this strong MSRO on the calculated 
SDR is now noticeable but still relatively small; the SDR 
is reduced compared to its MFA values by 12% for Fe 
and 22% for Ni. 



V. EFFECT OF THE LOCAL MOMENT 
REDUCTION 

Reduction of the local moment is a universal feature 
of itinerant magnets as revealed by spin fluctuations 
theories As discussed in Section Hill the local moment 
in Fe is very stable and changes only slightly in the para- 
magnetic state compared to zero temperature. Therefore, 
our calculations based on the ground-state value of the 
local moment agree well with experiment for Fe. Still, 
the SDR is sensitive to the local moment, and a small 
reduction of it noticeably improved the agreement with 
experiment at higher temperatures. Since the SDR was 
found to be insensitive to MSRO, it is reasonable to at- 
tribute the large overestimation of the high-temperature 
SDR in Ni to the neglect of the local moment reduction. 
Here we study this issue in detail. 

In the paramagnetic DLM state the local moment in 
Ni vanishes)^ but it is partially restored by longitudi- 
nal spin fluctuations. ^^'^^ Following the idea of separa- 
tion of low and high-energy fluctuations, we assume that 
the current-carrying quasiparticles near the Fermi level 
experience the averaged exchange-correlation field gener- 
ated by fast longitudinal spin fluctuations, and that this 
"mean field" is adequately represented by noncoUinear 
DFT with disordered local moments constrained to their 
square-averaged values. The atomic potentials are there- 
fore obtained using the fixed spin method^ with the 
value of the constrained local moment treated as an ad- 
justable parameter, which has a physical meaning and 
can be measured experimentally. Other approximations 
are, in principle, possible; for example, the longitudinal 
spin fluctuations can be explicitly included in the same 
noncoUinear DFT approach, i. e. they can be considered 
to be "slow" rather than "fast." Since the separation 
in slow and fast degrees of freedom is not well defined, 
we did not attempt to study the role of these additional 
fluctuations. 

The calculated paramagnetic SDR of Ni as a function 
of the local moment is shown in Table IHII As seen, SDR 
is very sensitive to the value of the local moment. Com- 
parison with experimental SDR shows that our predicted 
value of the square-averaged local moment in paramag- 
netic state of Ni is equal to 0.35/zb (using the more accu- 
rate basis set with Imax = 3). Unfortunately, we are not 
aware of experimental measurements suitable for com- 



TABLE III: Spin-disorder resistivity in fiQ-cm for paramag- 
netic Ni as a function of the fixed local moment. 2x2 super- 
cells and basis sets with Imax = 2 and Imax ~ 3 were used. 
Standard deviations of SDR due to limited disorder sampling 
are included. 
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"This value corresponds to unconstrained local moment of 0.63/jb . 



parison with this prediction. 

VI. DISCUSSION AND CONCLUSIONS 

Numerous previous studies^'^'^'^iii^ based on the s-d 
model concluded that SDR in the paramagnetic state is 
essentially proportional to J^j^S{S + 1) where S is the 
spin of the partially filled 3d shell. This dependence 
is easy to understand based on the Fermi golden rule 
with averaging over the initial states of the 3d spin. In 
our treatment based on noncoUinear DFT, the exchange- 
correlation field with randomized directions on different 
sites plays the role of the s-d Hamiltonian. However, con- 
trary to the s-d model, the 3d spin is treated classically, 
i. e. S is just a classical vector and not an operator. The 
Fermi golden rule in our case would give a paramagnetic 
SDR proportional to J'^^S^. Thus, if the S{S ± 1) factor 
were correct, noncoUinear DFT calculations would under- 
estimate the paramagnetic SDR by a factor {S ± 1)/S'. 
This factor is close to 2 for Fe and more than 3 for Ni. In 
reality, the calculated SDR agrees well with experiment 
for Fe and is overestimated for Ni (if the local moment 
reduction is not included). We believe that these results 
provide clear evidence against the S{S + 1) factor which 
appears if the local moments are treated as local atomic 
spins. Instead, the classical description of the local mag- 
netic fluctuations in the spirit of the DLM approach is 
supported by our results. We suggest that the itinerancy 
of the 3d electrons is crucial for this behavior. Qualita- 
tively, one can argue that the low-energy fluctuations in 
Fe or Ni on the scale of kT (which the resistivity is most 
sensitive to) are similar to classical rotations of the local 
moments rather than quantum fluctuations of localized 
spins. It would be interesting to investigate this issue for 
magnets with a varying degree of localization, including 
rare-earth systems. 

Some poorly controlled assumptions are involved in the 
extraction of Pmag from the experimental datai^ First, it 
is assumed that Pmag is constant in the wide temperature 
range above Tc where the total resistivity is linear in T. 
This assumption implies that the local moments (or at 
least their mean-squared average) are constant in this 
range. Spin fluctuation theories for itinerant metals show 
that the local moments may change with temperature 
above Te i^'^'^^i^^ Such change will contribute to the slope 
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of p above Tc, and hence the separation of Pmag from the 
phonon contribution would be inaccurate. 

On the other hand, it has been argued that the phonon 
contribution to the resistivity may be sensitive to spin 
disorder, because the latter may change the character of 
states at the Fermi level.—i^ In particular, in Ni the filled 
majority-spin d states may be lifted up to the Fermi level 
by spin disorder, thereby facilitating interband s-d scat- 
tering by phonons. This effect may therefore introduce 
an unusual temperature dependence of the phonon con- 
tribution, which makes spin disorder and phonon effects 
non-additive, even if the scattering rates themselves obey 
Matthiessen's rule. Since we have not studied this effect 
here, our comparison of SDR with experiment for Ni is 
incomplete. However, the phonon contribution can be 
expected to follow the Bloch-Griineisen form above 
with the electron-phonon scattering renormalized by spin 
disorder; therefore, the influence of spin disorder on the 
phonon contribution should not invalidate the procedure 
used for subtracting this contribution above Tc- 

In conclusion, we have calculated the spin-disorder re- 
sistivity of Fe and Ni in the whole temperature range 
up to Tc using both the mean-field approximation and 
the nearest-neighbor Heisenberg model to represent the 
canonical ensemble of classical spin configurations. We 
found that SDR is insensitive to the magnetic short-range 
order (MSRO) in Fe and Ni. The SDR in Fe depends lin- 
early on AP (T) which implies that the main effect of spin 
disorder is to introduce scattering, which is proportional 
to the variance of the random potential. For Ni the cal- 



culated temperature dependence is more complicated; at 
elevated temperatures close to Tc the SDR grows faster 
than expected. This faster increase of SDR may be ex- 
plained by the reduction of the exchange splitting which 
lifts the heavy bands up to the Fermi level, thereby in- 
creasing the scattering rate. The results for Fe are in very 
good agreement with experiment if the atomic potentials 
are taken from zero temperature and frozen, but for Ni 
the SDR calculated in this way is strongly overestimated. 
This disagreement is attributed to the reduction of the 
local magnetic moment in Ni. Comparison with experi- 
mental SDR leads to a value of 0.35/iB above Tc, which 
may be compared with experiment. 
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